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Geometrical methods for formation flying design based on the analytical solution to 
Hill's equations have been previously developed and used to specify desired relative motions 
in near circular orbits. By generating relationships between the vehicles that are intuitive, 
these approaches offer valuable insight into the relative motion and allow for the rapid 
design of satellite configurations to achieve mission specific requirements, such as vehicle 
separation at perigee or apogee, minimum separation, or a specific geometrical shape. 
Furthermore, the results obtained using geometrical approaches can be used to better 
constrain numerical optimization methods; allowing those methods to converge to optimal 
satellite configurations faster. This paper presents a set of geometrical relationships for 
formations in eccentric orbits, where Hill's equations are not valid, and shows how these 
relationships can be used to investigate formation designs and how they evolve with time. 


I. Introduction 

A great deal of attention in recent years has been focused on various ways to improve our ability to describe the 
motion of one vehicle about another (see Alfriend and Yan * 1 ). Broucke 2 , Alfriend and Ginr’. Lawden 6 , Melton 8 , 
and Tschauner and Hempel 1 1 all derived equations of relative motion valid for vehicles in nearby eccentric orbits and 
Alfriend and Gim’ and Schweighart and Sedwick 10 included first-order perturbation effects in their relative motion 
equations. This paper focuses on formation design in unperturbed, eccentric orbits using geometrical relationships 
derived from the relative motion equations described in Ref. 2, which are completely explicit in time. Note the 
relative motion equations developed in Refs. 3, 6, 8, and 1 1 all use true anomaly as the independent variable or are 
in the form of a series expansion in eccentricity. 

Previous geometrical methods for formation design and station-keeping, such as those proposed by Lovell and 
Tragesser 7 and Sabol et al 9 , are based on the analytical solution to Hill’s equations 4 , which limits their applicability 
to near-circular reference orbits. Recently, Hughes 5 published a numerical optimization method for formation 
design in eccentric orbits that maximizes a mission performance metric. While numerical optimization methods are 
necessary to maximize the performance of any mission, if used exclusively, they have several limitations. First, 
numerical optimization methods require a realistic initial guess (or guesses) of the optimal satellite configuration to 
avoid numerical instability problems. Developing a realistic initial guess requires some knowledge of the relative 
motion and the feasibility of a desired satellite configuration. Second, it is difficult to characterize how the 
properties of the relative motion are related to desired formation parameters. This relationship is typically explored 
in a brute-force manner: the relative motion properties for many different values of a design parameter are examined 
until the relationship between the two is understood. 

Alternatively, geometrical methods can be used to quickly design satellite configurations that achieve specific 
mission requirements. They offer an intuitive understanding of the relationship between the properties of the 
relative motion and desired formation parameters and the ability to visualize the relative motion without a 
simulation, a valuable tool early in the design process. Thus, the two methods compliment each other if used in 
conjunction. A geometrical design method that is not limited to near circular reference orbits is, therefore, desired. 
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Beginning with the recently published analytical solutions to the linearized equations of relative motion in an 
eccentric reference orbit provided in Ref. 2, we have developed a set of geometrical relationships for formations in 
eccentric orbits. These relationships allow for the use of intuitive design and station-keeping methods, similar to 
those discussed in Refs. 7 and 9, to be used to specify formation parameters. As an added benefit, these 
relationships directly relate initial conditions to Keplerian elements, eliminating some of the coordinate 
transformations necessary in other geometrical methods. 

The remainder of this paper develops the set of geometrical relationships for formations in eccentric orbits, 
shows how these relationships can be used for formation design, and, finally, examines the validity and accuracy of 
the proposed design methods for three test formations. 

II. Linearized Equations of Relative Motion in Elliptic Reference Orbit 

Recently, Ref. 2 presented a time-explicit solution to the linearized equations of motion of nearby objects in 
eccentric orbits. The linearization assumes only that the two objects are proximate in position and velocity, but does 
not require any assumptions about the eccentricity of the orbit. The equations are described in terms of the relative 
position and velocity of a deputy (Ref. 2 refers to this as the chaser) with respect to a chief (called the target in Ref. 
2) expressed in a curvilinear Local-Vertical-Local-Horizontal (LVLH) coordinate frame. 

This frame moves with the chief and will be referred to as the RAC (radial, along-track, and cross-track) frame 
in this paper. Figure 1 shows the system. The R axis is aligned with the position vector of the chief and the C axis 
is normal to the orbital plane of the chief. The A axis is perpendicular to the other two axes and its direction is 
defined by C x R. The R axis is a linear dimension; motions along the A and C axes are measured on a sphere that is 
locally tangent to the orbit of the chief at every instant in time (Fig. 1). The radius of the sphere is equal to the 
instantaneous position of the chief. The rotation rate of the frame is time varying. The use of this frame, instead of 
a more conventional Cartesian frame, dramatically improves the accuracy of the following analysis. 

The relative equations of motion in terms of the radial position (pc) and velocity, along-track position (y) and 
velocity, and cross-track position (z) and velocity, are given as 


x-29 y-9 2 x-9 y 
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The time varying coefficients in Eqs. (1) - (3) are expressed as 

2/fesin v 


where v is the true anomaly of the chief. 

In Ref. 2, the solution to Eqs. (1) and (2), the coplanar motion of the deputy, is derived by referencing the 
relative position and velocity to relative Keplerian elements. Specifically, Ref. 2 considers the semimajor axis a, 
eccentricity e, argument of perigee co, and mean anomaly at epoch Mo. Because we are also interested in out-of- 
plane motion, we consider the inclination i and the right ascension of the ascending node Q as well. Indicating the 
parameters of the orbit of the chief with no subscripts, we have 

O deputy — O + Aci, £ deputy — O + Ae, CO deputy — CO+ ACO, M deputed — M q + AM , 
i deputy ~ t'TAl, D deputy — £2 + AD, 
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where Aa , Ae , Aco , AM, A/, and AQ represent small deviations from the Keplerian elements of the chief. The exact 
solution to Eqs. (1) - (3) is expressed in terms of the Keplerian element differences as follows 
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z = r sin(<y + v)Az' — r sin i cos(<y + v / )Af2. 


( 5 ) 

(6) 


Note in Ref. 2 the expression for the along-track position, y, does not include a dependence on AO by definition 
because only coplanar motion of the deputy is considered and, for the out-of-plane motion to be zero, AO (and A;) 
must be zero. To understand what causes the along-track dependence on AO, consider the case when AO y 0 so that 
O, d e p U ty 7 = O, as illustrated in Fig. 2. In this case, the argument of perigee of the chief and deputy are measured in 
slightly different planes and this causes an offset between the chief and deputy in the along-track direction. 
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III. Geometry of Relative Motion in Elliptical Reference Orbit 

Beginning with Eqs. (4) - (6), we have established a set of geometrical relationships that describe, in an 
accessible way, the relative motion in an eccentric orbit. The goal was to create a simple set of descriptors, such as 
those used by Vallado 12 , which describe the motion of a deputy about a chief in an eccentric orbit. 

Because we are primarily interested in establishing stable formations, we first eliminate the possibility of secular 
growth in the separation between vehicles by constraining the energy of both orbits to be equal (i.e, A a = 0). We are 
left with 



ae sin v 

x = -a cos vAe + AM, 

( 7 ) 

y = 

f \ 2 

o + sinv , Ae+ — y/l-e 2 AM + r(A<y+cos/AD), 
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z = rsin(<y + i/)A; — rsinicos(<y+ v)A£l 

( 9 ) 


Note that all the terms in Eqs. (7) - (9) have an implicit dependence on time. With some algebraic manipulation, we 
can write Eqs. (7) - (9) as 

x = Csin(v / -^ 0 ), (10) 

y = c cos(v - y/ 0 ) ~ D cos(e + Y 0 ) ■ + y cm , 0 1) 

z = Gsm(E+(p 0 )+z cm , (12) 
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where E is the eccentric anomaly of the chief, y cm is the along-track center of motion, z cm is the cross-track center of 
motion, and the constants are defined as 
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The geometrical structure of Eqs. (10) - (12) can be used to investigate formation designs and how they evolve with 
time. 


IV. Formation Design 

Inspection of Eqs. (10) - (12) reveals that, in general, the relative motion in an eccentric orbit will not be an 
ellipse in any plane. This is a significant departure (and complication) from the analytical solution to Hill’s 
equations 4 , where the motion in the radial/ along-track plane is always a 2x1 ellipse with the major axis in the along- 
track direction. We will show later that elliptical motion is possible in the along-track/cross-track plane, but only 
when certain constraints are imposed on the motion (elliptical motion is also possible in the radial/along-track 
plane). The non-elliptical motion is due to the eccentricity of the reference orbit. In an eccentric orbit, E ^ v, except 
at perigee and apogee, and, thus, the terms that vary with eccentric anomaly have a different time dependence than 
those that vary with true anomaly. 

Five design parameters must be specified in Eqs. (10) - (12) in addition to the constraint that Aa = 0 to eliminate 
the possibility of secular growth between the vehicles. Three of the parameters, Ae, Aco, and AM, only affect the 
coplanar motion (i.e., radial/along-track motion) and A i only affects the out-of-plane motion (i.e., cross-track 
motion). Both the coplanar and out-of-plane motion is affected by AO. 

In the following sections, three types of formations are examined: (1) an along-track formation, (2) a follower 
formation, and (3) an along-track/cross-track formation. The design of each formation is described in detail below. 

A. Along-track Formations 

Purely along-track motion is achieved by requiring the radial and cross-track positions always be zero. Thus, C, 
G, and z cm must be zero in Eqs. (10) - (12). This is achieved by requiring Ae = AM = Ai = AO = 0. The only 
parameter to specify in an along-track formation design is Aco, which allows us to specify a desired along-track 
separation at a particular time in the orbit. For example, we can create a desired separation of p 0 at perigee by 
setting 


A (0 A f ~ 


Po 

a( 1-e) 
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or at apogee by setting 


t±(0 A F - 


Pa 

a( 1 + e) ’ 


where p A is the desired separation at apogee. 

The minimum separation, p+ n , between two vehicles in an along-track formation occurs at perigee ( V^ n = 0- 
deg) and is equal to 


Pmin =aAO) 4F (\-e). 


(14) 


Similarly, the maximum separation, p^ x , occurs at apogee ( = 1 80-deg) and is equal to 


pZl = aAa AF {l + e). 


(15) 


B. Follower Formations 

A follower formation is a formation of vehicles that share the same groundtrack. For two vehicles to occupy the 
same groundtrack, they must pass over the same set of points on the Earth. The important consideration here is that 
the nodal separation, AO, must account for AM and the rotation rate of the Earth, W e . In particular. 


„ AM 

AO pf = -W PP- 


(16) 


where n is the mean motion of the reference orbit (a detailed derivation of this relationship is provided in Ref. 9). It 
should be noted that Eq. (16) is only applicable when A a = 0. If A a P 0, then n for each vehicle is different and Eq. 
(16) is no longer valid. 

All the vehicles in a follower formation must have the same eccentricity, inclination, and argument of perigee if 
they are to share the same groundtrack. Thus Ae = A i = Aco = 0. The only free parameter in a follower formation 
design is AM, which represents a desired separation at a particular time in the orbit. We can create a separation of p 0 
at perigee by setting 


AM FF 



1 + e ( l-e)cosiW e 

2 

cos 2 &>sin 2 i W 2 (l - e) 2 

yyjl-e 2 n j 
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(17) 


where the ± indicates whether the second vehicle is ahead (+) or behind (— ) the reference vehicle. Substituting Eq. 
(17) back into Eq. (16) gives 


-p KPo 
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(18) 


C. Along-track/Cross-track Formations 

A formation whose motion is entirely in the along-track/cross-track plane is achieved by requiring the radial 
position always be zero. Thus Ae = AM = 0 and we are left with three design parameters to specify: Aco, Ai, and 
AO. 
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All along-track/cross-track formations resemble one of two types of motion in the along-track/cross-track plane. 
If the phase of the cross-track motion, tp 0 , equals 90-deg, then Eqs. (10) - (12) describe a straight line with slope, s, 
equal to 


\s 


a 

e(cos iAQ + Am) 


The sign of s is determined by the position of the vehicle at perigee; if the vehicle is in the bottom half of the plane, 
s is positive and if the vehicle is in the top half of the plane, s is negative. 

If (p 0 does not equal 90-deg, then Eqs. (10) - (12) describe an ellipse in the along-track/cross-track plane centered 
at [o(cos /AO + Am), -aeoc\. 

To simplify our analysis, we restrict the vehicles to linear motion, and, thus, from Eq. (12) we have the 
constraint 


Pacf = cos mAi ACF +sin/ sin mAQ. ACF = 0. (19) 

The remaining two design parameters characterize the size and shape of the formation at a particular time in the 
orbit. If we define the along-track and cross-track position at perigee to be y 0 and z 0 , then we can immediately 
identify the other two conditions. From Eq. (10) we have, 


cos /AO ACF 


+ A ® ACF ~ 


To 

a ( 1 - <?) 


and, similarly, from Eq. ( 1 1 ) we have 


sin mAi ACF 


-sin/cosftAO^ CF 


z o 

a( 1 - e) 


(20) 


(21) 


Special care must be taken when solving Eqs. (19) - (21) in equatorial (/ = 0 or 180-deg) and polar (/ = 90-deg) 
reference orbits, as well as when perigee of the reference orbit is located at a nodal crossing (co = 0 or 180-deg). We 
solve one such case here. If the reference orbit is inclined (/ ^ 0 or 180-deg) and perigee is located at a nodal 
crossing (co = 0 or 180-deg), then Eqs. (19) - (21) are solved by 


ACl 


ACF 


a(l-e)sin; 


Am 


ACF 


y 0 sin / + z 0 cos / 
a(l - e)sin / 

A * acf = 0. 


(22) 

(23) 

(24) 


The minimum separation between any two vehicles in an along-track/cross-track formation is a function of 
the phasing angle between the vehicles. The phase, of one vehicle with respect to another is 


r f 


= tan 


f 


\yoj 
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If z 0 /yo is less than or equal to e 1/2 , then the minimum separation between the vehicles, p ^ F , occurs at perigee 
( v m\n = 0-deg) and is equal to 


P 


ACF 

min 



+ Z n 


Otherwise, p^f F occurs at 


V ACF 
y min 


-1 


= COS 


( eyl ^ 


v z o j 


and is equal to 


p ACF 

A^min 


y o z o 


(e + l) 




2 2,2 

e y o + z 0 


(25) 


(26) 


(27) 


The maximum separation between any two vehicles, p AF [ , always occurs at apogee ( = 1 80-deg) and is equal 
to 


P 


ACF 

max 


1 + e 
1-e 


l 


2 , 2 
y o + z o • 


(28) 


Inspection of Eq. (27) reveals one important condition for along-track/cross-track formations: if y 0 = 0 then 
p m j,, = 0 and a collision will occur at V min . Therefore, vehicles in along-track/cross-track formations cannot be 
aligned in the along-track direction at perigee. 

V. Test Formations 

This section examines the accuracy of the formation design methods discussed in the preceding sections. An 
example of each formation type is considered: (1) an along-track formation, (2) a follower formation, and (3) an 
along-track/cross-track formation. The motion of a deputy about a chief is examined for each test formation. The 
size of each test formation at perigee (i.e., the separation between the chief and deputy) is 1-km. Because the size of 
the relative motion is proportional to the size of the formation at perigee, the relative motion for a different sized 
formation is obtained by simply scaling the results in this section by the desired size of the formation at perigee. 

Some of the characteristics of the relative motion are lost when the nonlinear equations of motion are linearized 
in the derivation of Eqs. (10) - (12). To evaluate the accuracy of Eqs. (10) - (12), a set of Keplerian elements is 
computed for each deputy using the design methods presented in this paper. The Keplerian elements of each deputy 
are converted to Earth-Centered-Inertial (ECI) coordinates, which are then propagated separately using simple two- 
body dynamics. These ECI coordinates serve as the ‘true’ position and velocity of the deputy. The true RAC 
position and velocity of the deputy are computed from the ECI coordinates using the following conversion, 


true 'deputy 
y true — true ’ 

z — yQ z 

true ' true ’ 

^ true ~ 'deputy ^ ’ 

V — rO y + Y0 y 

y true ' w true ' true’ 

Ztrue = rd tme + r ^r„e> 
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where x true , y trU e, and z tme are the true radial, along-track, and cross-track position of the deputy and 
* true’ y true* an d z true are the true radial, along-track, and cross-track velocity of the deputy. The details of the ECI to 
RAC conversion are summarized in Appendix A. 

The position linearization error, e p , is defined as the norm of the difference between the true RAC position and 
Eqs. (10) - (12), 


S„ = 


yl (Xtrue ~ X f + Gw ~ )’Y + ( Z true ~ Z Y 


(29) 


Similarly, the velocity linearization error, e n is defined as the norm of the difference between the true RAC velocity 
and the derivatives of Eqs. (10) - (12), which are found in Appendix B, 

£ v = VOw ~xf + (ytrue - yf + Gw - if ■ (30) 

Equations (29) and (30) are a measure of the information lost in Eqs. (10) - (12). It might be tempting to assume 
that Eqs. (29) and (30) occur independently of one another, but they are coupled. Both are a function of how nearby 
the deputy is to the chief in position and velocity; the key simplifying assumption in the derivation of Eqs. (10) - 
(12). This assumption is the main source of the linearization error and allows us to make an important observation: 
if you increase the size of the formation by an order of magnitude, you increase the linearization error by two orders 
of magnitude. Thus, the linearization error is proportional to the size of the formation squared. 

The orbit of the chief in all three test formations is the reference orbit. Table 1 lists the Keplerian elements of 
the chief. The chief appears stationary in each formation because it is always at the origin of the RAC frame. 


Table 1. Keplerian elements of chief in test formations. 


Parameter 

Chief 

Perigee radius ( r p ) 

16,072-km 

Apogee radius (r fl ) 

68,120-km 

Semimajor axis (a) 

42,096-km 

Eccentricity (e) 

0.6182 

Inclination (i) 

10-deg 

Right ascension of ascending node (O) 

0-deg 

Argument of perigee (oj) 

0-deg 

Mean anomaly at epoch (M 0 ) 

0-deg 


A. Along-track Test Formation 

The deputy only has an along-track component of motion in this formation. Using Eq. (13) and the values listed 
in Table 1, A coaf is 0.00357-deg. All the other Keplerian element differences are zero. 

The evolution of the formation in the radial/along-track plane is shown in Figure 3. The oscillatory motion of 
the deputy is illustrated in this figure. Figure 4 shows the true along-track position and velocity of the deputy. The 
true radial and cross-track position and velocity are zero within machine precision. In Fig. 4, the minimum 
separation between the vehicles is 1-krn and the maximum separation is approximately 4.24-km. Equations (14) and 
(15) predict the true minimum and maximum separation within machine precision. 

The position and velocity linearization errors for the deputy are always zero within machine precision. This is 
the result of the simple, linear motion of the deputy and the use of the curvilinear coordinate frame. 

B. Follower Test Formation 

Both vehicles share the same groundtrack in this formation, with the deputy following the chief. Using Eqs. (17) 
and (18) and the values listed in Table 1, A M FF is -0.00081-deg and AQ FF is 0.00081-deg. The other Keplerian 
element differences are zero. 
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The groundtrack of both vehicles is compared in Figure 5. The difference between the latitude and longitude of 
the chief at time t and the latitude and longitude of the deputy at time t + AM/77 is zero within machine precision and, 
thus, the chief and deputy share identical groundtracks. 

Figure 6 shows the true radial, along-track, and cross-track position and velocity of the deputy. Most of the 
motion of the deputy is contained in the radial/along-track plane, with only a slight cross-track motion. This is due 
to the small inclination of the orbit of the chief. As the inclination of the orbit of the chief increases, so does the 
magnitude of the cross-track motion. 

Figure 7 shows the position and velocity linearization errors for the deputy for 5 orbital periods. In Fig. 7, both 
the position and velocity linearization errors are periodic and bounded. The maximum position linearization error is 
approximately 0.022-m, which is less than 0.006-percent of the minimum separation between the vehicles. The 
maximum velocity linearization error is approximately 0.014-mm/s. In Figure 8, the size of the formation at perigee 
is varied from 1-krn to 100-krn. The position linearization error is not significant for most applications, never 
exceeding 1 -percent of the minimum vehicle separation for formation sizes less than 100-km. 

C. Along-track/Cross-track Test Formation 

The motion of the deputy in this formation is linear and contained in the along-track/cross-track plane. Figure 9 
plots AQ acf an d A co A cf, which were computed using Eqs. (22) and (23) and the values listed in Table 1, as a 
tunction of the phase of the deputy. The other Keplerian elements differences are zero. 

The motion and evolution of the deputy in the along-track/cross-track plane for several different phases is shown 
in Figures 10 and 11. The straight-line motion of the deputy is illustrated in both figures. Figure 12 shows the true 
along-track and cross-track position and velocity of the deputy as a tunction of the phase of the deputy. The radial 
position and velocity is zero within machine precision. As we can see in Figs. 10 - 12, purely along-track motion 
results when the phase of the deputy is 0 or 180-deg (i.e., when z 0 = 0) and purely cross-track motion results when 
the phase of the deputy is 90 or 270-deg (i.e., wh cn y 0 = 0). Purely along-track motion is discussed in section A and 
the vehicles collide wheny 0 = 0, so the reminder of this section does not discuss the 0, 90, 180, or 270-deg cases. 
They are included on figures whenever necessary as envelopes of the relative motion. 

Figure 13 shows the true separation between the vehicles as a tunction of the phase of the deputy. The 
maximum separation is independent of the phase and is approximately 4.24-km. The error in Eq. (28) is small, 
typically less than 10" 4 -percent of the true value. If the phase of the deputy is less than or equal to tan _1 (e 1/2 ), then 
the true minimum separation is approximately 1-km. The minimum separation predicted by Eq. (25) is also 
generally within 10" 4 -percent of the true value. If the phase of the deputy is greater than tan"'(e 12 ), then the true 
anomaly at which the minimum separation occurs predicted by Eq. (26) and the minimum separation predicted by 
Eq. (27) are both typically within 0.01 -percent of the true value. 

The position and velocity linearization errors for the deputy as a function of the phase of the deputy for 5 orbital 
periods are shown in Figure 14. The interesting result is deputies with supplementary phases, such as 60 and 120- 
deg, have different linearization errors. At first, this seems inconsistent: the magnitude of the position and velocity 
for deputies with supplementary phases appears to be the same in Fig. 12 and, therefore, both deputies should have 
the same linearization error. 

However, consider for a moment two deputies, one with a phasing of 60-deg (labeled Dl) and one with as 
phasing of 120-deg (labeled D2). Figure 15 shows the cross-track components of the position and velocity 
linearization errors for Dl and D2. Examination of Figs. 14 and 15 reveals that the cross-track components account 
for almost all the position and velocity linearization errors. This is due to a lack of symmetry in the cross-track 
motion of Dl and D2 that is not captured by Eqs. (10) - (12). According to Eqs. (10) - (12), the magnitude of the 
cross-track position and velocity of Dl and D2 are equal at every instant in time. This, however, is not realistic 
because Dl is ahead of D2 in its orbit. This causes Dl and D2 to have slightly different cross-track positions and 
velocities at every instant in time. Thus, the position and velocity linearization errors for Dl and D2, and all 
deputies with supplementary phases, are different. This effect is also present in the along-track direction, but is 
more noticeable in the cross-track direction because of a greater lack of symmetry between the two orbits in this 
dimension. The along-track components of the linearization error are generally an order of magnitude smaller then 
the cross-track components. 
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The position and velocity linearization error for both deputies is again periodic and bounded in Fig. 14. The 
maximum position linearization error is always less than 0.4-m, which is approximately 0.15-percent of the 
minimum separation between the vehicles. The maximum velocity linearization error is always less than 0.08- 
mm/s. The size of the formation at perigee is varied from 1-km to 100-km in Figure 16. The position linearization 
error exceeds 1 -percent of the minimum separation at a formation size of approximately 7-km. 
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Figure 3. Evolution of deputy in along-track formation in radial/along-track plane. 



Figure 4. True RAC coordinates for along-track formation. 
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GroundTrack 



Figure 5. Groundtrack for follower formation. 



Figure 6. True RAC coordinates for follower formation. 
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Figure 7. Linearization error for follower formation. 



Figure 8. Maximum linearization error for follower formation. 
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Keplerian Element Differences for Along-track/Cross-track Formation 



Figure 9. Keplerian element differences for along-track/cross-track formation. 



Figure 10. Motion of deputy in along-track/cross-track formation in along-track/cross-track plane. 
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Figure 11. Evolution of deputy in along-track/cross-track formation in along-track/cross-track plane. 


True RAC Coordinates for Deputy in Along-track/Cross-track Formation 
Position 






Figure 12. True RAC coordinates for along-track/cross-track formation. 
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Figure 13. Inter-vehicle separation for along-track/cross-track formation. 
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Figure 15. Cross-track components of linearization error for along-track/cross-track formation. 




Size of Formation at Perigee (km) 

Figure 16. Maximum linearization error for along-track/cross-track formation. 
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VI. Conclusions 

This paper developed and evaluated a set of geometrical relationships for formation flying design in eccentric 
orbits. Geometrical design methods are valuable for two reasons: (1) they provide a rapid method to design 
formation geometries and constraints to meet specific mission requirements and (2) they offer valuable insight into 
the characteristics of the relative motion. The geometrical relationships presented in this paper are based on the 
analytical solution to the linearized equations of relative motion in an eccentric reference orbit discussed in Ref. 2. 
Previous geometrical methods for formation design have been limited to near circular reference orbits, but we have 
avoided this restriction. 

Three types of formations were analyzed: (1) an along-track formation, (2) a follower formation, and (3) an 
along-track/cross-track formation. Our formation design method proved to be a useful tool not only in the design of 
appropriate formation constraints, but also in understanding the characteristics and properties of the relative motion. 
The position linearization error was less than 0.4-m and the velocity linearization error was less than 0.08-mm/s for 
a 1-km formation at perigee. 

The next steps will be to (1) evaluate the influence of perturbations and higher-order gravity terms on the 
motion, (2) use the results presented in this paper as a starting point for numerical optimization methods, and (3) 
employ Eqs. (10) - (12) in a relative navigation filter. 
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APPENDIX A: Rotating from ECI to RAC 

If the ECI coordinates of the chief and deputy are known such that 

r= \ x ECI y ECI z ECl\ *‘=\x E ci ECI Z ECI ] > 

^ deputy — EC! V ECI z ECI ] deputy 5 ^ deputy ~ ECI ECI z ECI ] deputy’ 5 

then the RAC position and velocity of the deputy with respect to the chief is compute using the following 
conversion, 


r ^ r r • r _ 
R= — , r =tz\ r r 

M i —i 3 


Irl’ 


c = r~zi’ c = 0 

r Xr 


A=CxR, A = CxR , 


R 


deputy 


' deputy 


^ deputy 'deputy 'deputy _ 


_ ’ deputy _ 


'deputy 'deputy 

^ deputy 


3 ' deputy ’ 


0 V = sin 


-\ f - - \ • R deputy ‘ A + R deputy ' ^ 

deputy Oy ~ ^ ’ 

COS 6',. 


8 _ = sin 


-i th r\ R - Rde P u( y ' C 

V s - deputy W’ u z 

cos 0 7 


X = \ r deputy]- \r\, X = 


r deputy ' r deputy V • V 


'deputy 


y= rO v , y=-rr-d v + rd , 


z= rd,, z='-^-8 z + r 8,, 


where x, y, and z are the radial, along-track and cross-track position of the deputy with respect to the chief, 
and x, y, and z are the radial, along-track and cross-track velocity of the deputy with respect to the chief. 
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APPENDIX B: The velocity equations 

The velocity equations are obtained by taken the time derivative of Eqs. (7) - (9), 


k = n sin v-\/l — e 2 
n-Jl — 


'a* 


Ae + ewcosv 




i AM, 
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1+-^ 
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r 2 
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.... aen sin v , 

I AM + — A a>+ 




an 




t (sin i (sin 6 + e sin + (cos 6 + e cos co)ki). 
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